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Abstract. Filiz et al. (2008) proposed a model for the pattern of defaults 
seen among a group of firms at the end of a given time period. The ingredients 
in the model are a graph G = (V, E), where the vertices V correspond to the 
firms and the edges E describe the network of interdependencies between the 
firms, a parameter for each vertex that captures the individual propensity 
of that firm to default, and a parameter for each edge that captures the joint 
propensity of the two connected firms to default. The correlated default model 
can be re-rewritten as a standard Ising model on the graph by identifying the 
set of defaulting firms in the default model with the set of sites in the Ising 
model for which the {±l}-valued spin is +1. We ask whether there is a suitable 
continuous time Markov chain (Xt)t>o taking values in the subsets of V such 
that Xq = 0, X r C X a for r < s (that is, once a firm defaults it stays in 
default), the distribution of Xx for some fixed time T is the one given by the 
default model, and the distribution of Xt for other times t is described by 
a probability distribution in the same family as the default model. In terms 
of the equivalent Ising model, this corresponds to asking if it is possible to 
begin at time with a configuration in which every spin is —1 and then flip 
spins one at a time from —1 to +1 according to Markovian dynamics so that 
the configuration of spins at each time is described by some Ising model and 
at time T the configuration is distributed according to the prescribed Ising 
model. We show for three simple but financially natural special cases that this 
is not possible outside of the trivial case where there is complete independence 
between the firms. 



1. Introduction 

Investors are exposed to credit risk due to the possibility that one or more coun- 
terparties in a financial agreement will default] that is, not honor their obligations 
to make certain payments. Some examples of default are a consumer or business 
not making a due payment on a loan, a manufacturer or retailer not paying for 
goods already received from a supplier, a bond issuer not making coupon or prin- 
cipal payments, or an insolvent financial institution not returning deposited funds 
to its customers upon demand. 

Some credit risk is present in virtually any financial agreement, and a key ingredi- 
ent in its satisfactory management is a model that produces a sufficiently accurate 
probability for a given default event. Consequently, there is a large theoretical 
and applied li terature on this topic |BR02I IBOW02I IDSMI ldSR.041 IGief)4bl ISch04l 
ZP07, Wag08 . Roughly speaking, models of default lie on a spectrum between the 
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structural and reduced form ones. For the example of a firm defaulting on its debt 
obligations, a structural model might include explicit descriptions of the dynamics 
of the firm's assets, capital holdings and debt structure, whereas a reduced form 
model would not seek to incorporate the details of the actual mechanism by which 
the firm is led to default but rather it might typically be something of a "black box" 
that treats the time of default as a random time with an associated exogenous in- 
tensity process having a rather simple structure characterized by a small number of 
parameters which may have little direct economic interpretation. Although struc- 
tural models are perhaps theoretically more satisfying because in principle they 
provide a means of testing how well the factors that cause default are understood, 
they are often perceived as being too complex and parameter-rich for them to be 
fitted adequately: defaults are uncommon and even firms within the same sector of 
the economy can be quite heterogeneous, so there can be insufficient "independent 
replication" upon which to base statistically sound parameter estimates. 

The difficulty of modeling default probabilities is compounded for complex fi- 
nancial instruments such as collateralized debt obligations (CDOs) and other struc- 
tured asset-backed securities that are constructed by, in essence, bundling together a 
group of borrowers. It is then no longer sufficient to determine the default probabil- 
ity for a single "firm" - rather, it becomes necessary to model the joint probabilities 
that various subsets of firms in the basket will default, and it is usually not ap- 
propriate to treat the defaults of different firms as statistically independent events. 
The most obvious reason for this absence of independence is that firms are subject 
to the same background economic environment. Moreover, situations such as the 
interconnectedness between manufacturers, their parts suppliers and the retailers 
who sell their products can cause problems for one firm in such a network to spread 
to others via a process that is usually described as contagion. A small sampling of 
the substantial empirical and modeling work on this phenomenon of correlated de- 
fault is [DLOlbl IDEoTal IHW01I ISSOll IZhoOll IFMMI iGleMl ILG031 IGie04al IGW04I 
IGW06I IDDKS071 IEGG071 iJZlTFl lYuTffl lESMl iGGDMl ICDGHlOj . 

In this paper we investigate a particularly appealing class of models for correlated 
default in |FGMS08j (see also |MV05[ IKMFI06] . where special cases of this model 
were introduced). The basic model in |FGMS08] does not attempt to describe the 
time course of defaults for some group of firms. Rather, it is a one period model 
that gives the probability any given subset of the firms will have defaulted at some 
time during a prescribed time interval. 

The ingredients of the model in FGMS08 are a finite (undirected, simple) graph 
G with vertex set V and edge set E and two vectors of parameters a = (a v ) v< zv G 
M. v and (3 = ((3 e )eeE £ ■ Each vertex v G V represents a firm and the graph 
structure provided by the edges is intended to capture the network of interdepen- 
dencies between the firms. Write I v , v G V, for the indicator random variable of 
the event that firm v defaults; that is, /„ takes the value 1 if firm v defaults and 
the value otherwise. The probability of a given pattern e = (e v ) v ^v € {0, 1} V of 
defaults is 



(1.1) 



P{J, 



e v , v G V} 



-exp(F(e)), 



where the Hamiltonian H is given by 



(1.2) 




w 



uev 



{v,w}£E 
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and the partition function Z is the normalizing constant that ensures the sum over 
{0,1}^ of the probabilities is one. The parameter a u , u G V, is clearly some 
measure of the individual propensity of firm u to default. The parameter /3{„ jl0 }, 
{v, w} G E, captures in some way the dependence between the defaults of firm 
v and firm w: if this parameter is positive, then the joint default of both firms 
is favored, whereas it is discouraged when the parameter is negative. We write 
V{G, a, /3) for the distribution of the random binary vector /. 

Note that if we set Y v = 2I V — 1, v G V, then (Y v ) veV G {±1} V and for 
{o~v)vev G {±1} V we have 

P{Y V = a v , v eV} = P{I V = (1 + a v )/2, v G V} 

= ^ ex p I \ X) a «( 1 +°'tt) + \ / 3 {«,«'}( 1 + cr «)( i + cr «>) 

V uGV {v,w}eE 

= -~ exp j 7„cr„ + ^ 5{ V:W ya v a w 

\u£V {v,w}EE 

for suitable parameters (j v )vev G ^ V and (8 e ) e ^E S K s and a corresponding 
normalization constant Z . Thus, the random vector of spins {Y v ) v ^v is described 
by the usual Ising model associated with the graph G — (V, E). 

It is shown in FGMS08 that this class of correlated default models is as flexible 
as one could possibly hope: if J is an arbitrary {0, l} y -valued random variable, then 
there is a choice of the parameters (a v ) v£ v ancl {Pe)e£E such that I u has the same 
distribution as J u for all u g V and for all {v, w} G E the pair (I V ,I W ) has the same 
distribution as (J V ,J W ). Moreover, it is observed in FGMS08 that it is possible 
to fit such a model to data using existing techniques such as iterative proportional 
fitting, various convex optimization techniques, or a number of other "off-the-shelf" 
numerical optimization methods suitable for large-scale computation. 

A significant drawback of the class of models in |FGMS08] is that they don't 
provide a description of the time dynamics of default: they just give the probability 
that a given subset of firms have defaulted during some fixed time period without 
saying anything about the distribution of the times at which the defaults occurred. 
If we let [0,T] be the time period of interest, then we would like there to be a 
{0, l} y -valued stochastic process (I(t))o<t<T such that 

• I v (t) = 1 if and only if firm v € V has defaulted by time t, so that 
ii,(0) = and the sample paths of (I v (t))o<t<T are right-continuous and 
non-decreasing (once a firm defaults it does not "undefault" ) , 

• #{i> G V : I v (t) ^ I v (t—)} < 1 for any t G [0,T] (two or more firms do not 
default simultaneously), 

• I(T) is has distribution V{G, a, /3). 

Furthermore, since V(G, a, ft) is supposed to be an appropriate description for the 
pattern of defaults during [0, T], it is reasonable to require that 

• I(t) has distribution V(G,a(t), (3(t)) for suitable parameters a(t) and /3(i) 
when < t < T. 

In this paper we investigate whether such a process exists within the simplest 
and perhaps most natural class of models, namely the time-homogeneous Markov 
chains. Re-cast in the language of the equivalent Ising model, we are thus asking if 
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it is possible to begin at time with a configuration in which every spin is —1 and 
then flip spins one at a time from —1 to +1 according to Markovian dynamics so 
that the configuration of spins at time T is distributed according to a prescribed 
Ising model and at all other times the configuration is described by some Ising 
model. 

We can certainly construct such a chain if f$ = 0, so that V(G, a, (3) = V(G, a, 0) 
is the distribution of avector (I v ) ve y of independent {0, l}-valued random variables 
with 

F{I V =Q}= 1 

1 + exp(a„) 

We simply takes the processes (I v (t))t>o to be independent, with 

¥{I v (t) =0} = exp(-A„t), 
where the jump rate X v is chosen so that 

exp(-A t ,T) = — — — — -. 

1 + exp(a„) 

Thus, X v — ^log(l + exp(a„)) and I(t) has distribution V(G, a(t), 0), where 

exp(-A.„i) = exp ( — — log(l +exp(a 1 ,)) 



so that 



1 I I f4\\ I rri 

1 + exp(a„(t)) V T 

a v (t) = log ((1 + exp(a 1 ,))T - l) 



for < t < T. 

After establishing some general facts in Section [21 we investigate in Sections [31 
[H and [5] whether it is possible to construct a time-homogeneous Markov chain for 
non-zero j3 in the following cases: 

(I) G is the complete graph Kn in which there are N vertices with each vertex 
connected to every other one, a u — a v for u, v 6 V and /3 e = /3f for 
e,f&E; 

(II) G is the complete bipartite graph Km,n in which V is partitioned into two 
disjoint subsets V and V of cardinality M and N such that every vertex in 
V is connected to every vertex in V and there are no other edges, a u = a v 
for u, v G V and (3 e = /3f for e, / G E; 
(III) G is again the complete bipartite graph Km,n, ol u — a v for u,v G V, 
a u — a v for u, v e V and f3 e = j3f for e, / G E. 

In Model I there is complete symmetry: each firm has the same individual propen- 
sity to default and the interdependence between any two firms is the same as that 
between any other two. Model II and III both describe a situation in which there 
are two types of firms (say, for example, car manufacturers and auto parts suppli- 
ers) and there is only interdependence between firms of different types. In Model 
II all firms have the same individual propensity to default, whereas in Model III 
this propensity can depend on the type of the firm. 

We conclude in all three cases (with a minor technical restriction for Model 
III) that it is impossible to construct a time-homogeneous Markov chain with the 
desired properties unless /3 is zero; that is, unless the firms behave independently. 
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2. Generalities 

It will be notationally more convenient to identify a vector e = (e v )veV S {0, 1} 
with the subset A = {v G V : e v = 1} C V and regard V(G 7 a,f3) as a probability 
measure on subsets of V rather than {0, 1} V . If we extend the definition of f3{ u v } 
by declaring that /3( UjlI } = when {u,v} ^ E and write /3{ Mj „} more simply as fi UVl 
then our Hamiltonian, now thought of as function defined on subsets of V, is given 

by 

(2.1) H(A):=Y / ^u+ ]T Puv 

uEA {u,v}<ZA 

If we write P H for the probability measure V(G, a, /3), then 

(2.2) F H ({A}):=±exp(H(A)), 



where 



BCV 



We are interested in the existence of a time-homogeneous Markov chain X = 
(X t )t>Q that has as its state-space the collection of subsets of V and has the fol- 
lowing properties, where we write Q(A, B) for the jump rate from state A to state 
B: 

• Q{A, B) — unless B = A U {v} for some v A; 

• when X(0) = 0, the distribution of X(T) is P H ; 

• there are parameter vectors a(t) and j3(t) for < t < T such that if we set 

H t (A):=J2*u(t)+ E <M*)> 

UEA {u,v}cA 

then X(t) has distribution F Ht when X(0) = 0. 

If such a Markov chain exists, we say that the default model admits time- 
homogeneous Markovian dynamics. 

Write A — > B if _B = ^4u{t>} for some u ^ A. The Kolmogorov forward equations 
for the chain X with initial state become 

(2.3) j¥ Ht (B)= £ P Ht (A)Q(A,B)+P fff (B)g(B,B), 

A— >B 

where, as usual, we put Q(-B, B) := — J2c^b Qi^i C)- 

Denoting the partition function associated with the Hamiltonian H t by Z t := 

EceB e ^ (C) ' we have 

— P fft (B) = 

dt at Z t 

_ Z t e H ^H[(B)-Z[e H ^ 



and thus (|2.3[) can be re- written as 

(2.4) Z t e H ^H' t (B) - Z' t e H ^ = ^ Q(A, B)e Ht ( A) Z t + Q(B, B)e Ht ^Z t 



A^B 
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To further simplify notation, set Rb = —Q(B,B). Because Ht{%) — 0, equation 
(|2~4|) for B = is simply 

—Z t = —Rt$Z t . 

We require 

1 = limP fft (0) = lim— , 
tj.0 tio Z t 

and so 

(2.5) Z t = e R ^. 
Substituting (1231) into (j2~l|) gives 

(2.6) H' t {B) = Q(A,B)e Ht ^- H ^ + R $ - R B , 
and hence 

ueB {u,v}CB 

(2-7) / \ 

= ^Q(B\ {«}, B) exp -a u {t) - ^ *M*) + ^ ~ 
«efi \ «efl\{o} / 

For ueV, set Q u = Q(0, {«}) and R u := i? {tl} = -Q({w}, {u}). Equation (I2T7)) 
for £? = {u} is 

(2.8) a^(i) = Q M e- Q "W + i? - 

Hence, by the method of variation of parameters (also called variation of constants), 

(2.9) att(t )= log 
and 

i?0 — R u 



(2.10) c/ u (t) = 



1 - e -(Rft-Ru)t 
when i?0 ^ R u . If i?0 = then 

(2.11) a u (t)= log (Q u t) 
and 

(2.12) a' u (t) = \. 

Note that each function a u , u G V, is completely determined by the rates Q u = 
Q(0jW), #0 = E^ev^'M)' and i?„ = J2vev\{u} Q(i u }> { u > v }), and hence 
the vector of functions (a„) u£ y is completely determined by the collection of rates 
{Q(0, {u}) : u e V} U {Q({u}, {w, u}) : u, v G V, u # «}. 
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For u,v G V, set Q uv := Q({u},{u,v}) and R uv := R{ u>v } = —Q({u, v}, {u, «}). 
Equation ()2.7[) for B = {u,v} is, upon substituting from (|2.10j) . 

- a4(t) - a4(t) - i? u „ + i? 

_ QvU R$ ~ Ru —/3 U v(t) I ^ Z ~Pvu{t) 



R$ — Ru 



i?0 — Ry 



1 — e -(.Re-Ru)t 1 _ e -(R»-Rv)t 
+ i?0 — i? ul , 

when i? ^ i?„ and R ^ R v . Analogous results hold when R = R u or R = R v . 
Recall that f3 uv (t) = /3r UtV \(t) — P vu (t), and so (|2.13l) is an ordinary differential 
equation for the function /3 UV if we treat the rates of the Markov chain as given. In 
particular, the two vectors of functions (a u )ueV an d (Pw)u,veV,v^v ar e completely 
determined by the collection of rates {Q(0, {u}) : u G V} U {Q({u}, {u, v}) : u, v G 
V, u 7^ v} U {Q({u, v}, {u, v, w}) : u,v,w G V, u ^ v ^ w ^ u}. 

In principle, we could attempt to find values for these rates such that (a u (T)) ue v 
and ((3uv(T))u,vev,u^v have the required value, substitute the resulting values of 
a u (t) and f5 uv (t) into ^7} (using (|?J0|l or (j2~T2"]) for the values of a' u (t) and (gUBD 
or its analogues when i? = R u or R — R v for the values of P' uv {t)) and hope to 
either find values for the remaining rates so that (|2.7|) holds for all B C V or show 
that this is impossible no matter what our initial choice of rates was. This seems to 
be a rather forbidding task in general, but we are able to carry it out in the three 
special cases described in the Introduction. 

3. Model I: complete symmetry 

Recall Model I from the Introduction. The graph G is Kn, the complete graph 
on N vertices for some N, and there are functions a and j3 such that 

| a u (t) = a(t) for all u G V 

I Puv(t) = P(t) for all u, v G Vu ± v. 



(3.1) 



Proposition 3.1. Model I with N > 4 admits time-homogeneous Markovian dy- 
namics if and only if the firms default independently. 

Proof. We observed in the Introduction that the general default model admits Mar- 
kovian dynamics when firms default independently. So we need to establish a con- 
verse for the special case of Model I with N > 4. 

Suppose that a collection of rates exists such that (|2.7|) holds for all subsets B. 
When #B > 1, (27} becomes 



(3.2) 



P'(t) 



J2Q(B\{u},B) 
ueB 
+ R$ — Rb- 



exp(-a(t)-(#S-l)/3(i)) 



If we average Q3.2j) over all ( ■ ) choices of sets B with #2? = k for some k > 1 and 
set 

Ra, 0<£<N, 



E 

ACV,#A= 



s 
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we get the equations 

ka'(t) + (§P'{t) = (Ac - X k ) 



(3.3) 



+ X k -i AT : e-^-f-VM, l<k<N. 
N — k + 1 



Note that X e > for < I < N - 1 and X N = 0. 
Equation (|3.3p for k = 1 and k — 2 yields 



(3.4) 



a'(t) = (A -A 1 ) + ^ e ^( t ) 

P'{t) = (A - A 2 ) + Ai^e-^W-^*) - 2(A - Ai) - 2^ e -*< t >. 



Substituting the values for a'(i) and /?'(£) from (J3T4J) into (|3.3p gives a system of 
equations of the form 

(3.5) Xk ^ e-^-iJ/sW-Ai— e'^ = a k e a ^ +b k , 1 < k < N, 

for appropriate constants a k and b k , 1 < k < N , that depend on the constants A^, 
< I < N. 

We claim that the continuous function /3 is constant. Suppose that this is not 
so. Note that a k can be non-zero for at most one value of k € {1, . . . , N}, because 
if a k , ^ and a k „ ^ for 1 < k! < k" < N, then 



\,, , fL e -(fc'-l)/3(t) _ V _2_™ _ b 



\,„ , hi e -{k"-l)P{t) _ x _2_ _ l 



A: 

and letting i vary over an open interval J such that the image {(3(t) : t G J} 
contains an open interval we would conclude that two polynomials of different 
degrees coincided over an open interval. Because N > 4, we thus must have a k = 
for some k > 3. Observe for such a k that 

x k -i — - — e-( k -ww _ Xlr J_ e -m = b 

N~k+1 N-l 

and again we would conclude that two polynomials of different degrees coincided 
over an open interval. Therefore, the function j3 must be a constant, say [3*. Of 
course, (3* is the pre-specified value for (3(T). 

We now show that (3* — 0. Equation (|3.3|l now becomes 



ka\t) = (A - A fe ) + A fc _ l jv _ fc fc+i e-"W-( fc -^ 



and hence, by (|3.4j) . 
(3.6) 



(Ao-A^ + ^e 



(A -A fc) + A fc _ l¥ -^— e -W-(*-^-. 



Each side of Q3.6P is first degree polynomial in e for every k £ {1, . . . , N}. It is 
apparent from the differential equation in (|3.4[) that the function a is not constant 
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(indeed, we solved this equation explicitly in (|2.9p ). Consequently, the coefficients 
of these two polynomial coincide and hence 

r fc(A - Ai) = (A - A fe ) 

(3-7) I fc^£-A t fe c -(k-m> 

for 1 < k < N. Re-arranging (J377J) , we conclude that 

A fc = fc(Ai - A ) + A 

for 1 < k < N and 

_ (N-k) fc/r 
*k — jz A e 

for < k < N — 1. Because N > 4, this is impossible unless j3* = 0. □ 
4. Model II: Two classes with common individual propensity to 

DEFAULT 

Recall Model II from the Introduction. The graph G is Km,n, the complete 
bipartite graph with vertex set the disjoint union V = V U V, where V has M 
vertices, V has N vertices, and there are functions a and (3 such that 

( a u (t) = a(t) for all u e V 
\ Puv(t) = fi{t) for all u € V, v € V. 



(4.1) 



Proposition 4.1. Model II with M > 3 or N > 3 admits time-homogeneous 
Markovian dynamics if and only if the firms default independently. 



Proof. As in the proof of Proposition 13. li it suffices from the remarks made in 
the Introduction about the general model to show that if the model admits time- 
homogeneous Markovian dynamics, then the firms default independently. 

Symmetry considerations similar to those in the proof of Proposition 13.11 show 
that if (|2.7[) holds for some choice of jump rates, then there are constants A^ n and 

Xln „, < m < M and < n < N, with A^„ = for < n < N, \l l N = for 
< to < M, and A^ n and n strictly positive otherwise such that 

(to + n)a'(t) + mn(3'(t) = r- A^ n - A^„ 

m \-> p - a (t)- n p(t) 

'Vn— l,n° 



( 4 - 2 ) M-m + 1 



n 



.. 1 'Vi,n-1 



\t -a(t)-m0(t) 



where we set r := A J + A^ and adopt the convention that A^j „ = 0, < n. < AT, 
d A^ _! = 0, < to < M. We leave t 
Setting (m,n) = (1,0) in (l4~2j) gives 



and A^ _ x = 0, < to < M. We leave the straightforward details to the reader 



(4-3) a '( t ) =r _(A£ + Al i0 ) + -M 

Similarly, setting (m, n) = (0, 1) in (|4.2p gives 

(4-4) r-(A^ + Aj, 1 ) + ^ 



A °^° e -"(t). 
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In particular, we have the identity 



An 



^0,0 



A 0,0 

M N ' 

Setting (m, n) = (1, 1) in (14.21) and substituting in the expression for a'(t) from 
(g31) gives 



(4.5) 



(r - (A^o + Al.o) + 




-a(t)-/3(t) 



-«(t)-/9(t) 



2a' (t) 



M 



Further substituting the above expressions for a'(t) and f3'{t) from (|4.3p and (I4.5[) 
into (|4.2[) for a general pair (to, n) leads to a system of equations of the form 



(4.6) 



o«(*) 



-n/3(t) 



< to < A/, < n < A 
where the various coefficients are given by 

t 



m/3(t) 



(to + n — 2mn)(r — ATq — A{ ) + TOn(r — A{ : — Aj\) 



At, „ + A 



(to + n — 2mn) ■ 
toA™_i ~ 



m.n 
An^r 



M 




Note that c m) „ > whenever to > and d m ^ n > whenever n > 0. 

We claim that the continuous function (3 is constant. Assume without loss of 
generality that M > 3 and suppose that the function j3 is not constant. 

Consider the two cases: 

(i) One of 02,1 or 031 is zero. 

(ii) Both 02,1 and 03,1 are non-zero. 

Case (i) is impossible, because for either to — 2 or to = 3 we would have 
= d m ^e~ m PW + (c rrii i — e m .i)e~P^ — b m .\ and conclude that either a quadratic 
or cubic polynomial was identically zero on some open interval. 

Turning to Case (ii), note first that if a m ^ n 7^ we have 



(4.7) 



c e-™^*) + d e - m0 M - e e~W - b 



and so we would have a quadratic and a cubic polynomial that agreed on an open 
interval. 



CORRELATED DEFAULT 



11 



Therefore, the function /? must be a constant, say f3* . Of course, /3* is the 
pre-specified value for f3(T). We now show that (3* = 0. 
Equation (|4.2[) becomes 



(4.8) 



M - m + 1 



a(t)-nP* I ro,n-l -gffl-mff* 

JV-n + 1 



For a fixed pair (m, n), each side of (|4.8|) is a first degree polynomial in a{t) and, 
since a(t) is continuous and non-constant, we can equate coefficients. If we also 
record our boundary conditions and conventions from above we arrive at the fol- 
lowing system of equations for < to < M and < n < N 



(4.9) 



(m + n) 
M 



A^n — 



M-m + 1 



-71/3* 



N-n + 1 

(to + n)(r - A^ - A^ ) = r - A^ n - A^ 
Al 



-A T p" 



m/3* 



A 



t 



m , — 1 
A-X.n = 0. 



Because A^J w = \ N = 0, we see from the second equation of (|4.9[) for (m, n) = 
(M, AT) that 



(4.10) 



r-Ka- aT , 



M + N' 



and we can substitute this value into the second equation of (|4.9I) for general (to, n) 
to conclude that 



+ >L,n = r ( 1 - |^ 1 n - .U. ()-_-,,-_- A 



Setting m = in the first equation of (|4.9[) and noting the fifth equation, we get 

N-n+1 



A 0,7l-1 



M 



-Xon, 1 < n < N. 



A similar argument leads to 



, ^ M — m + 1 , __. 
^-i.o = ~r X o,o, 1 < m < M. 



M 
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Combining these observations, we arrive at the following system of equations 



(4.11) 



n ^m,n-l (■ 

■ e K ' 

N — n + 1 

\m + n)e n0 ' 



mX, 



t 

m — l,n 



'0,0 

M 



M-m+1 

m(M + N-m-n + l) 



1 < m < M, < n < N, 
m + n 



M + N 



M-m + 1 



< m < M, < n < N, 



Xm-i = °- < m < M, 

X^ l n = 0, < n < N, 

A^„ = 0, < n < N, 

A!„ n = 0, < m < M, 



N — n 



M 

M - m 
M ' 



-Aq^o, < n < N, 



< m < M. 



Note that if A^„ and A^ „ satisfy (gUJ), then they also satisfy (14. lip . It will 
thus suffice to show that if /3* ^ 0, then there do not exist A^ „ and n satisfying 

Km . 

Setting (m,n) = (1,1) in the first equation of (|4. 1 1[) gives 
Alo Aq.i _ Ao!o f p _ M + N -I 



2e^ 
t 



M 



1 -° ~ M ,0 ' 



N M M 
while the fifth equation of (jlTl) forces A$ ^ = - 1). Thus, 

(4-12) Al )0 = ^A^ (2e^-l). 

If instead we set (m, n) = (2, 0) in the first equation from (14. lip , we obtain 
(4.13) ^ 

Comparing (14.121) and (|4.13p . we conclude that 

2e r -1 = 1, 

and hence j3* = 0. □ 

5. TWO CLASSES WITH DIFFERENT INDIVIDUAL PROPENSITY TO DEFAULT 

Recall Model III from the Introduction. As with Model II, the graph G is Km,n, 
the complete bipartite graph with vertex set the disjoint union V = V U V, where 
V has M vertices and V has N vertices. Now, however, there are functions a, a 
and (3 such that 

a u (t) = a(t) for all u G V 
(5.1) { a v (t) = a(t) for all v e V 

Puv{t) = P(t) for all u e V, v € V. 
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Proposition 5.1. Consider Model III with M > 4,JV > 3 or M > 3, N > 4. 

Suppose that the prescribed values of a(T) and a(T) are distinct. If the prescribed 
value of f3(T) is non-zero and sufficiently small, then the model does not admit 
time-homogeneous Markovian dynamics. 



Proof. Another symmetry argument similar to those in the proofs of Proposition l3.ll 
and Proposition 14.11 shows that if (12.71) holds for some choice of jump rates, then 
there are constants A^„ and A^ „, < to < M and < n < N, with A^ n = 

for < n < N, A^ N = for < m < M, and A^„ and Ajj ^ strictly positive 
otherwise such that 



ma'(t) + na'(t)+mn/3'(t) = r-(X^ tn + X"l ntn ) + 
(b - Z) + M-m+l 6 + N - n + i e 



where we set r := A J + A^q and adopt the convention that Az*i n = 0, < n < iV, 
and A^ _! = 0, < to < M. 

Applying (15.21) with (m,n) = (1,0) and (m,n) = (0, 1) gives 



(5.3) 



d'(i) = r - (A* + A-) + 



M 



-&{£) 



N 



Similarly, applying (|5.2p with (m,n) = (1,1) and then substituting in the expres- 
sions for a(t) and a(t) from (|5.3|) gives 



= r - (Aj 4 + A-) + ^e-«W-/>W + ^° e -*«-™ - tf(t) - a'(t) 
( 5 - 4 ) = -r - (Al iX + A-) - (Xl, + Aft) - (A^ + A^) + *£ e -*M-««) 

A l,0 affl-ffffl _ A oto -aft) _ Vo Q (t) 

M N 



Substituting the expressions for the a, a and from f)5.3[) and (J5T4J) into (|5.2 
for general (to, n) produces a system of equations of the form 



a e a(t) +fr +r P «(*)-«(*) + W 4. P P <*(*)-<i(*)-/3(*) 

( 5 - 5 > _ f p -nfi(t) , a(t)-&(t)-mff(t) 
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where 



: (m - mn)(r - X[ - A^ ) + (n - mn)(r - A^ - A^) 
+ mn(r - \\ A - A^) - r + A^„ + A^„ 
A^ n 

: (to — mn j ■ 



(n — mn) 



M 

A 0,0 

iV 



drn.n — 



n — mn- 



&m,n — Tim- 



M 
A i,o 



fm,n — 



mA 



m — l,n 



M - m + 1 



«AL 



AT - n + 1 

Observe that because A^ n is strictly positive for 1 < to < M — 1 and 1 < n < N 
and n is strictly positive for 1 < m < M and 1 < n < N — 1, both /,„.„ and 
<?m, n are strictly positive for 1 < to < M and 1 < n < iV. 

We claim that the continuous function /3 is constant. Assume without loss of 
generality that M > 4, N > 3 and suppose that the function /3 is not constant. 

Re-arrange (|5.5p to get 



(5.6) 



„ _|_P p-"(*)_Lp p-a(t) p -/9(0 _ „ p-a(t) p -™/3(t) 



„<*(*) 



— fm,n& 



-n/9(t) 



d e _,8(t) -6 



Because / m! n is strictly positive for 1 < to < M and 1 < n < AT, there is an open 
interval J such that the image : f £ J) contains an open interval and the 

right-hand side of (|5.6|) is non-zero for all t £ J, 1 < to < M and 2 < n < N, and 
hence the same is true for the left-hand side. 

Taking (|5.6[) with the indices (m,n) replaced by another pair (i, j), we see that 
if t £ J, 1 < m, i < M and 2<n,j <N, then 



, -n/3(t) _ J -/3(t) _ b 



(5.7) 



-&{t) 



-a(t)„-/3(t) 



fee 



-a(t) -i/3(t) 



V- , e -J/3(i) _ j. -/3(t) _ b . . 



Re-arranging (|5.7|) gives 
(5.8) p(e-^W; m, n, i, - q(e~ m ; m, n, i, j), 

where 

777* , Thji^jj . — (^ijj ^i.j ^ ^ 9i,j ^ ) \f 771,71% n% 

(Cm,n H - 6m,n^ Qm,n% )\fi,j^ di^jZ 
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and 



q(z;m,n,i,j) := a mtn (f itj z J - d itj z - b^) - a ij (/ m: „z r 



Suppose now that 2 < to, i < M and 2 < n,j < N. The leading term of the 
polynomial p(z;m,n,i,j) is -gi^fm.nZ^ 71 if i + n > to + j and g m ,nfi,jZ m+0 if 
i + n < m + j (recall that f m .n,gm.n, fi,j>9i,j are all strictly positive). Therefore, 
by taking a subinterval of J if necessary, when i + n ^ m + j we may suppose 
that J retains the properties required of it above and, moreover, that both sides of 
(|5.8p are non-zero for all t £ J. In particular, either a mjI1 7^ or a^j 7^ and the 
polynomial q(z;m,n,i,j) has degree either n or j when n j. 
Consider two 4-tuples (to' \n' ,i' , j') and (to", n", i", j") with 

2 < m',m",i',i" < M 



(5.9) 



O ^- / // ■/ ■// ^ AT 

2 < n , n , j , j <jv 

■'1 ' ' 1 ■' 

z + n /m + j 

i" + n" ^ m" + j" 



We conclude from (|5.8|) that 

(5.10) p(z; m', n', i',j')q(z; m" , n" , i", j") = p(z; m", n" , i",j")q{z; m', ri, j') 

for all z in an open interval. The left-hand side of (|5.10|) is a polynomial in z of 
degree either ((i'+n') V(m'+j'))+n" or ((z'+rt') V(m / +j'))+j", whereas the right- 
hand size has degree either ((z" + 7i")V(to"+j'")) + ?i' or ((i" + n") V(m"+/')) +/. 
For (m',n',i',f) = (2,2,2,3) and (m" ,n" ,i" = (4,2,4,3) we have 

4 ^ 5 = m' + j' 
m" + j" 



i + n 

i" + n" = 6 ^ 7 
n' = 2 ^ 3 = f 
n" = 2 ^ 3 = j" 
((z' + ra') V (to' - 



3')) 



= 5 



whereas the 



2 = 7 

{[i 1 + n') V (to' + j')) + j" = 5 + 3 = 8 

((i" + rz") V (to" + j")) + n' = 7 + 2 = 9 

((<" + n") V (to" + j")) + / = 7 + 3 = 10, 

and so the possible degrees of the left-hand side of (|5. 10[) are 7 and 
possible degrees of the right-hand side are 9 and 10. 

Therefore, the function f3 must be a constant, say f3* . Note that j3* is just the 
pre-specified value for /3(T). We now show that j3* = 0. 

For the moment, consider Model III with P(T) = 0, so that the function (3 must 
be identically zero and the firms evolve independently. In this special case, we know 
from the Introduction that 

(5.11) exp(d(<)) = (1 + exp(a(T)))T 

and 



1 



(5.12) 



exp(d(<)) = (1 + exp(d(T)))T - 1. 
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It follows from the linear independence of the functions exp(ci-), . . . , exp(c?r) when 
ci, . . . , Ch are distinct that in this case the functions exp(d), exp(d) and exp(d + a) 
are linearly independent when d(T) ^ d(T). 

Now return to the case of a general value for /3* = f3(T). Equation (|5.2p becomes 

ma'(t) + na'(t) = r- (A^„ + A^J 
(5.13) yy, \ y „\t 

, " LA m-l,n c -a(t)-nP' , /tA m,rt-l e _a(t)-m^* 

M - m + 1 iV - n + 1 

We can solve (|5.3[) for the d and d as in Section [2] to get 



(5.14) 



a(t) = log ( — *f° - (e^UX _ i) 
\M(r -A{ i0 -A^ ) , 

fl(t) = log f (eC-^-Vx)* - 1)) 

\N(r — Aj i — A^ x ) J 



Substituting ([534]) into ([535)) gives 

m - A lo - A l?o) + M r - A 0J - A o!l) + ( A m,n + A m,n) ~ r 

^ f _ m^. 1|B e-^' \ M(r-\j -\^ ) 
M M-m + l ) w r P (r-\tn-\? )t 



(5.15) V M M-m+l j A^ (e^- A i."-W - T 

m,n— 1 



'nAL nAl,„_ ie -^*\ AT(r -A^-A^V) 



AT 



AT-n + 1 y xl^-^-^ -1) 



= 0. 



It follows from the observations above and a compactness argument that if 
a(T) ^ d(T) and (3* = f3(T) is sufficiently close to zero, then the functions exp(d), 
exp(a) and exp(d + d) are linearly independent. Suppose that this is the case. 
Equation (|5.15[) is of the form 

a + be- &W + ce- &(t) = 

for suitable constants a, b, c, and hence a = b = c = 0. Thus, 

mA qo _ fflA m-i,/"" r = Q 
M M-m+l 

nA 0,0 ' tA m,n-l e 





N N -n+ 1 

m ( r - A l?0 - A Lo) + n ( r - Kl - A S : l) + ( A m,» + X ™,n) ~ T = 

for (to, n) € {0, . . . , M} x {0, . . . , N}, and so, after some algebra, 

A ?o^V r + Xlo^e^' =r-m(r- X{ ) n(r X^ A^) 

for (m, n) G {0, . . . , M} x {0, . . . , AT}. In particular, considering (to, n) = (fc, fc) for 
< k < M A JV leads to a system of the form 

Ake kf3 ' + Be kf3 ' + Ck + D = 

for suitable constants A,B,C,D with A > and i? > 0. A straight line can 
intersect the graph of the function t n> Ate 13 1 + Be 13 * at most twice if (i* > and 
at most three times if (3* < 0, and since M A JV > 3 we must have (3* = 0. □ 
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